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We discuss one possible method for suppression of the nonlinear effect of a parametric oscillatory 
instability in a Fabry-Perot cavity of the next generation of gravitational wave detectors by deter¬ 
mination of optimal radii of curvature of interferometers’ mirrors. Varying the radii of curvature of 
the mirrors we can change the frequencies of the optical modes of the Fabry-Perot resonators and 
distributions of optical fields in optical modes, which, in turn, can change the number of unstable 
modes, leading to parametric oscillatory instability. The determination of the optimal values of the 
radii of curvature of mirrors to reduce the number of unstable modes by thermal lensing is fulfilled. 


INTRODUCTION 


Currently, the European Gravitational Observatory is 
developing a project of a third generation gravitational 
wave interferometer that would allow solving more am¬ 
bitious problems than the gravitational detectors such as 
LIGO, VIRGO, GEO and others. The future third gen¬ 
eration instrumentation for searching for gravitational 
waves, which was named the Einstein Telescope(ET), will 
have a sensitivity several orders of magnitude higher than 
that of second generation gravitational wave detectors. 
The details of this project are still in the development 
stage and the first data on the Fabrv-Perot cavities that 
form the ET arms were outlined in [l|i3( . 

Distance L between the Fabry-Perot cavity mirrors will 
also be increased to 10km compared to 4km in LIGO, 
mirror radius will be increased to 31cm, which is about 
a factor of 4.25 larger than the laser-beam radius on the 
mirrors, thus making optical mode losses due to diffrac¬ 
tion very small. To reduce the level of exposure to seismic 
influence, it is planned to place the third generation grav¬ 
itational wave detector in a horizontal tunnel with two 
10-km arms underground(see all other parameters of ET 
in Table G}. 

In full scale terrestrial interferometric gravitational 
wave antennae like ET the undesirable effect of paramet¬ 
ric oscillatory instability (PI) in the Fabry-Perot (FP) 
cavities may take placeQ. Due to the high optical power 
of up to IT ~ 3MW circulating in the arm cavities a 
transfer of optical energy driving mechanical oscillations 
of the mirrors is likely. This behaviour might disturb the 
locking performance of the interferometer and can cause 
a substantial decrease of the antennae sensitivity. 

The effect of PI describes the excitation of the mirror’s 
elastic modes fed by the optical energy in the resonator. 
Thus, the optical main mode (TEMoo) with frequency 
wo effectively drives the elastic motion of the mirror with 
frequency w m as well as a second optical mode with fre- 
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quency wi which is called Stokes mode. Especially at the 
resonance wo—wi = w m the described energy transfer be¬ 
comes very efficient. The coupling between these three 
modes arises both due to the ponderomotive pressure of 
light in the main and Stokes modes and the paramet¬ 
ric action of mechanical oscillation on the optical modes. 
Above a critical value of optical power W c in the res¬ 
onator the amplitude of the mirror’s mechanical oscilla¬ 
tion rises exponentially as well as the optical power in 
the Stokes mode. 

In 2002 it was shown that the mirror’s oscillation could 
be damped by the same mechanisms. So the mechani¬ 
cal oscillation of the mirror could drive an optical anti- 
Stokes mode with frequency w a = wo + w™ resulting in 
the damping of the motion of the mirror [5!j. This will 
of course lessen the negative consequences of PI for an 
interferometer. But in [b] it has been shown that the 
consideration of anti-Stokes modes could not completely 
suppress the effect of parametric oscillatory instability. 
This calculation was based on the geometry of the LIGO 
interferometer including the power recycling mirror. Also 
the detailed analysis of Pis for the LIGO interferometer 
including signal recycling has been provided in HQ. 

The condition for the appearence of a parametric os¬ 
cillatory instability for the FP cavity may be written in 
a simple relation for the parametric gain 1Z Q: 
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Here c is the speed of light, L is the length of the FP 
cavity, m is the mirror’s mass, 7 m and 71 are the re¬ 
laxation rates of the elastic and Stokes modes respec¬ 
tively, W is the power circulating inside the cavity and 
A = wo — Wi — w m is the detuning value. Aj characterises 
the overlapping factor between the elastic and the two 
optical modes (main mode and Stokes mode). Aq and 
Ai represent the optical field distributions at the reflect¬ 
ing surface of the mirror for the optical main mode and 
the Stokes mode, respectively. f dr± corresponds to the 
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integration over this reflecting surface and f dV - over 
the mirror volume V. The vector u describes the spatial 
displacement of the mirror’s elastic mode and u z is the 
z-component of u oriented along the cylindrical axis. 

As we can see from Eq. 0 the parametric instability is 
a threshold effect and it takes place if the optical power 
W in the main mode of the FP cavity is bigger than 
a threshold value W c . The value W c depends on the 
detuning A between the optical and elastic modes and 
the overlapping factor Ai. In particular, W c takes its 
minimum value if the detuning is less than the relaxation 
rate of the Stokes mode A <C 71 and the overlapping 
factor is large enough. 

Analyzing Eq. 0 the authors in [Bj pointed out that 
the presence of the anti-Stokes mode may depress or even 
exclude Pis. Indeed, including the effect of an anti-Stokes 
mode with frequency u}\ a results in a negative gain 
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where 71 a is the relaxation rate of the anti-Stokes mode 
and A a = wi a — ujq—< jj m is the detuning value for the anti- 
Stokes mode. This damping effect has to be considered 
in a PI search. The effective parametric gain for each 
elastic mode has the following form lZ tot = TZ + TZ a - 

In particular, we need detailed information about pairs 
of Stokes modes and clastic modes which may be possi¬ 
ble candidates for parametric instability. It is known that 
Stokes and anti-Stokes modes may be analytically calcu¬ 
lated for Gaussian beams. In contrast, elastic modes can 
only be calculated numerically. 

Imperfections of the mirrors from the cylindrical shape 
(such as flats and suspension ’ears’) cause both splitting 
of originally mechanical degenerate modes into doublets 
and an additional frequency shifts of the elastic inodes. 
If, in turn, the difference between doublet frequencies 
is large enough, i.e. greater than the relaxation rate of 
the Stokes mode new Pis could emerge. Thus, the ap¬ 
pearance of doublets may increase the possibility of PI 

In order to predict parametric unstable combinations 
of Stokes and elastic modes we use an additional az¬ 
imuthal condition. Let us assume TEMoo as the optical 
main mode. If the elastic mode and the Stokes mode 
show an azimuthal dependence of ~ e l1 ^ and ~ e im< ^, 
respectively, the overlapping factor will vanish unless 
m = l. This behaviour results from the azimuthal in¬ 
tegration in Eq. 0 for Ai. It is worth noting that this is 
correct only if the cylinder center and the laser spot cen¬ 
ter coincide. For Laguerre-Gaussian optical mode LG33 
which uses as a pump in ET the azimuthal condition is 
\m±l\ = 3. 

In this Letter we analyse PI in ET at different val¬ 
ues of radii of curvature and find optimal ones for ef¬ 
fective suppression of PI. In section I we introduce all 


details of our PI calculations. In section II we deduce 
the dependence of the number of unstable modes on ra¬ 
dius of curvature(ROC) for the end test mass(ETM) and 
find optimal range of its ROCs. The ROC of input test 
mass(ITM) is fixed as in Table Q] 


I. PI IN ET 


The most important requirement for a PI search con¬ 
sists in the knowledge of the mechanical resonant fre¬ 
quencies and their elastic mode shapes. For elastically 
isotropic and cylindrical substrates there exist analytical 
solutions specialised to some restrictive conditions. One 
example is the well-known Chree-Lamb modes HM 1 
that are only applicable for certain ratios of the cylin¬ 
der’s radius and height and show no azimuthal defor¬ 
mation. Recent investigations using the superposition 
method allow to calculate the resonant frequencies and 
mode shapes of all axially symmetric modes of a cylin¬ 
der with high accuracy. In their work Meleshko et al. 
0 derived partial solutions of the equation of motion 
that satisfy the boundary condition of the free sample. 
But the only method to obtain the complete set of elastic 
mode shapes and frequencies for our PI calculation con¬ 
sists in a finite element modelling. For this purpose the 
FE package COMSOL® was used. In this Letter we pro¬ 
pose accuracy estimates to be about 1 % due to analytical 
solutions of Chree-Lamb-modes. 

For the mechanical relaxation rate 7 m we use the model 
of structural loss and express the relaxation rate with 
respect to mechanical loss (f>. Within this approximation 
we find a frequency dependent mechanical relaxation rate 

7 m ~ . (3) 

Note that this approximation is only valid for small me¬ 
chanical loss values. In the further calculation we use a 
constant value of (f> = 1 x 10~ 8 for fused silica at room 
temperature. 

Optical resonant frequencies and the spatial distribu¬ 
tion of their electric field on the surface of the mirror de¬ 
fine the possibility of PI. The spectrum of optical mode 
frequencies depends on the geometry of the FP cavity, 
which is realised in the interferometer arms of ET. 

The optical mode frequencies in a Fabry-Perot cavity 
can be calculated due to Kogelnik [2C| as 


WNnm = £ {N 0 TT + (2n + TO + 1) X 0 O ) , 
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There Nq = 0, ±1, ±2,... is the longitudinal index, n = 
0 , 1 , 2 ,... is the radial index and m = 0 , 1 , 2 ,... is the 
azimuthal index. R\ and i ?2 are the radii of curvature for 
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the ITM and ETM of the Fabry-Perot cavity(see details 
in Table |T|. The optical main mode circulating inside the 
cavity in ET is identified by n = 3 and to = 3(LG 33 
mode). 

Consequently, the transition from the optical main 
mode wo to another optical mode wi is accompanied by 
a frequency change of 

Q 

Aw = wi — wo = y [(2n + to) x <f>o — ttN] . (5) 
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Note that only negative frequency changes Aw < 0 can 
cause Pis as they belong to an optical Stokes mode. In 
contrast positive frequency changes belong to anti-Stokes 
modes, which can effectively suppress Pis and will be 
considered in a separate calculation. 

It is worth noting that not only the frequency detun¬ 
ing A but also the loss of optical energy influences the 
parametric gain. The latter is summarised in the op¬ 
tical relaxation rate 71 . On the one hand, energy loss 
occurs due to transmission of the mirrors and, on the 
other hand, due to clipping losses as a result of diffrac¬ 
tion . For our calculation we only consider the coefficient 
of energy transmission for the input test mass Titm as 
well as clipping energy losses on both mirrors T c u p . Both 
these coefficients create effective transmission coefficient 


TtotUM- 


This calculation shows an increase of optical loss for 
increasing mode orders due to the clipping loss fraction. 
Thus, the effect of parametric instability of high order 
optical modes becomes more and more unlikely (see 71 in 
Eq. (H])). In this work we only considered optical modes 
whose transversal orders fulfill the condition 2n+TO < 10. 


By the same time the mirror is excited by a PI includ¬ 
ing a Stokes mode it could as well be damped by the 
effect of PI due to feeding of an anti-Stokes mode. Sub¬ 
sequently the gain for the mirror oscillation could fall 
below the threshold of TZ = 1 and prevent parametric 
instability to occur for the given oscillation. 

In a detailed analysis we took the elastic mode of each 
candidate and calculated the parametric gain 7Z a of the 
anti-Stokes effect for optical modes satisfying 2 n + m < 
10. Finally we added all anti-Stokes gains lZ a and sub¬ 
tracted this value from the original parametric gain value 
TZ. This procedure lead to a total of nine parametric un¬ 
stable combinations. 

11 1 1 1 1 1 1 1 1 
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FIG. 1: The number of unstable elastic modes with IZtot > 1 
while tuning the radius of curvature of ETM 


II. RESULTS 

The calculation of parametric gain TZ follows from 
Eqs. (P), © and ©. As the elastic frequencies are or¬ 
ders of magnitude smaller than the optical ones we can 
use the relation wi = wo — w m ~ wo- Together with 
the dispersion relation of light and the wavelength of the 
optical main mode A 

wi/c ~ wo/c = 27r/A ( 6 ) 

we can simplify the expression for the parametric gain to 
_ 16-7T W 1 Ai 
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Using a MATLAB® code we evaluated this expression 
for mechanical modes up to a resonant elastic frequency 
of 30 kHz and for optical Stokes modes confined to mode 
orders of 2n + to < 10. With these numbers we obtain a 
total of nine possible PI candidates in ET for the men¬ 
tioned range of optical and mechanical modes. 


One possible method to minimize the number of unsta¬ 
ble modes is the changing of mirrors’ radii of curvature 
by local heating using additional laserfthermal lensing). 
This idea was firstly pointed out in @,0, 221 for LIGO in¬ 
terferometer. The thermal tuning of optical Fabry-Perot 
cavities has been proposed to reduce the parametric gain 
of parametric instabilities[22j]. The authors investigated 
the performance achievable for such tuning obtained by 
thermal actuation of the mirrors of the arm cavities. In 
our Letter we present the similar method for reduction of 
Pis and propose to vary the radius of curvature of ETM 
in one Fabry-Perot cavity. The range of these variations 
is from 5623 to 5663 meters. We took this interval into 
account because of making small changes in ET design. 

In FigCQthe dependence of number of unstable modes 
on the radius of curvature of ETM is plotted. In our cal¬ 
culations we assume that each non-axisynmietric elastic 
mode has sine and cosine component in azimuthal angle 
with practically the same frequencies. Therefore, each 
component gives an equal contribution to the paramet¬ 
ric gain IZtot- In Fig|T]we consider each unstable elastic 
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mode as one unstable mode having these two compo¬ 
nents. It is worth noting that the planed value of radius 
of curvature for ETM(I ? 2 = 5643m) is not optimal one 
because of large number of unstable modes (nine unsta¬ 
ble modes only in frequency range of elastic modes up to 
30kHz). Varying the radius of curvature of ETM R 2l for 
example, by thermal lensing(the R\ is a constant one) we 
can suppress the PI approximately by a factor of 2 (five 
unstable modes) if the radius of curvature of ETM mirror 
becomes slightly bigger. The prefered values of possible 
radii of curvatures lie in two ranges from 5648.3 to 5650.8 
meters and from 5651.7 to 5655.5 meters. Because of 
some inaccuracy in elastic parameters and numeric cal¬ 
culations the better interval for suppression of PI is a 
second long interval. In Table HJ all unstable combina¬ 
tions of elastic and optical modes with appropriate para¬ 
metric gains are shown(see Appendix A). If we do not 
take into account the influence of modes splitting some 
elastic modes will be stable. For example, for the ROC 
of ETM R 2 = 5623m elastic modes with the frequencies 
17091Hz and 25339Hz become stable(for R 2 = 5633m - 
20285Hz; for R 2 = 5643m - 13572Hz, 13701Hz, 25339Hz; 
for 5653m - 23387Hz; for 5663m - 28398Hz correspond¬ 
ingly). 

CONCLUSION 

In this Letter we have analyzed the possibility of the 
reduction of the unstable modes in ET by varying, for 
example, the radius of curvature of ETM mirror. In this 
case the frequencies of the optical modes of the Fabry- 
Perot resonators and distribution of their optical fields on 
the mirrors’ surfaces change, which, in turn, can change 
the number of unstable modes leading to parametric os¬ 
cillatory instability. 

Another important problem related to PI is the in¬ 
crease in the elastic modes’ density for high frequencies 
which is caused by deviations of the shape of real test 
masses of the ET from an ideally cylindrical shape(they 
have ’’ears” for suspension and symmetrically arranged 
flat cuts on the sides of the cylinders [Hi, |T(|). Such 
changes of the mirror shape result in a shift of the fre¬ 
quency spectrum of elastic modes and splitting of non- 
axisymmetric modes into doublets which increase the 
density of elastic modes by nearly twofold. These fac¬ 
tors, in turn, stimulate high probability of the onset of 
the PI in ET and change the conditions during which the 
PI is observed. 

It is worth noting that the mean interval between 
neighboring elastic modes near the elastic frequency 
u) m = 27r x 30000Hz can be estimated as Auj m ~ 
100s _1 [1, 0|. But in real detectors the mirrors have 
flats and suspension ’’ears” which increase the number 
of mechanical modes and, in turn, the probability of PI. 
A value of splitting gap between sine and cosine com¬ 


ponents for non-axisymmetric elastic modes can archive 
Su) m ~ 40s -1 . Neglecting the clipping losses(for optical 
modes with 2n + m < 10 these losses are negligible) we 
find relaxation rate of these modes to be 71 cs 55s -1 . 
Therefore, for elastic modes near the elastic frequency 
Lo m = 27r x 30000Hz (taking into account the splitting) 
the mean interval between neighboring elastic modes is 
approximately equal to the relaxation rates of the op¬ 
tical modes Au> m /2 — 71 . It means that elastic spec¬ 
trum becomes continious and practically for each elastic 
mode with frequency higher than 30kHz there exists the 
Stokes mode with small detuning A < 71 . For this elas¬ 
tic modes range the number of unstable modes will be 
practically the same for different radii of curvature. But 
only accurate calculation of overlapping factors is needed 
for searching of unstable elastic modes. 

Our results have stochastical character because of vari¬ 
ations in elastic parameters (Youngs modulus, density 
and others). Therefore, the values of radii of curvature 
can change and influence the number of possible Pis. On 
the other hand, the presence of power recycling cavity 
causes parametric gains IZtot to be hugely amplified[ 6 j. 
In this case the experimental results are needed. 


Acknowledgement 

The author kindly acknowledges the support of the 
Russian Foundation for Basic Research(grants nos. 11- 
02-00383-a and 14-02-00399-a). 


A. Parameters and detailed results of the 
calculations 

In this Appendix we present the numerical values of the 
parameters of ET that influence our calculations. As the 
substrates are made from fused silica we used material 
parameters according to this material at room tempera¬ 
ture. The numerical values are presented in Table H) The 
parametric gains for some values of radius of curvature 
for ETM are shown in Table [TT1 
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R 2 [m] 

Ri [m] 

fm [Hz] 

l 

n 

m 

TZtot 

5623 

5643 

8486 

0 

2 

3 

370.17 

5623 

5643 

8764 

0 

2 

3 

4.6 

5623 

5643 

13701 

1 

1 

2 

50 

5623 

5643 

13701 

1 

0 

4 

6.66 

5623 

5643 

16913 

2 

2 

1 

123.21 

5623 

5643 

17091 

2 

2 

1 

1.17 

5623 

5643 

20194 

1 

1 

4 

32.6 

5623 

5643 

20194 

1 

2 

2 

7.52 

5623 

5643 

20285 

1 

2 

2 

2.12 

5623 

5643 

25339 

2 

1 

1 

1.3 

5623 

5643 

25394 

2 

1 

1 

3 

5623 

5643 

28682 

1 

0 

4 

154.26 

5623 

5643 

28682 

1 

1 

2 

1.15 

5633 

5643 

8486 

0 

2 

3 

70.61 

5633 

5643 

8764 

0 

2 

3 

3.93 

5633 

5643 

13701 

1 

1 

2 

4.39 

5633 

5643 

16913 

2 

2 

1 

292.1 

5633 

5643 

20194 

1 

1 

4 

13.55 

5633 

5643 

20194 

1 

2 

2 

3.13 

5623 

5643 

20285 

1 

2 

2 

1.01 

5633 

5643 

23368 

0 

2 

3 

1.1 

5633 

5643 

25339 

2 

1 

1 

81.81 

5623 

5643 

28682 

1 

0 

4 

3.46 

5643 

5643 

8486 

0 

2 

3 

28.05 

5643 

5643 

8764 

0 

2 

3 

3.4 

5643 

5643 

12096 

1 

3 

2 

1.7 

5643 

5643 

13572 

1 

1 

2 

1.5 

5643 

5643 

16913 

2 

2 

1 

26.97 

5643 

5643 

20194 

1 

1 

4 

2.31 

5643 

5643 

21967 

3 

1 

0 

4.25 

5643 

5643 

23368 

0 

2 

3 

2.22 

5643 

5643 

25339 

2 

1 

1 

1.62 

5653 

5643 

8486 

0 

2 

3 

14.92 

5653 

5643 

8764 

0 

2 

3 

2.98 

5653 

5643 

16913 

2 

2 

1 

9.21 

5653 

5643 

23368 

0 

2 

3 

6.46 

5653 

5643 

23387 

10 

0 

7 

1.67 

5663 

5643 

8486 

0 

2 

3 

9.24 

5663 

5643 

8764 

0 

2 

3 

2.63 

5663 

5643 

16913 

2 

2 

1 

4.6 

5663 

5643 

23368 

0 

2 

3 

48.63 

5663 

5643 

23387 

10 

0 

7 

2.24 

5663 

5643 

28411 

3 

2 

0 

5.43 

5663 

5643 

28398 

1 

1 

2 

1.22 


TABLE II: Candidates for parametric instabilities. Given 
are the mirrors’ radii of curvature of Fabry-Perot cavity, the 
elastic resonant frequency f m = w m /( 2n), azimuthal num¬ 
ber of elastic mode (Z), the radial (n) and azimuthal (m) in- 
deces of the optical mode and the effective parametric gain 
IZtot = 1Z + 7 Z a - The contribution of possible anti Stokes 
modes is considered in the calculations. 








